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Abstract : 



We study wave equations with energy dependent potentials. Simple analytical models are 
found useful to illustrate difficulties encountered with the calculation and interpretation of 
observables. A formal analysis shows under which conditions such equations can be handled 
as evolution equation of quantum theory with an energy dependent potential. Once these 
conditions are met, such theory can be transformed into ordinary quantum theory. 
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1 Introduction 



Wave equations with energy dependent potentials have been known for long time. They ap- 
peared along with the Klein-Gordon equation for a particle in an external electromagnetic 
field [1] . In non-relativistic quantum mechanics, they arise from momentum dependent interac- 
tions, as shown by Green [2]. The Pauh-Schrodinger equation represents another example [3, 4]. 
A Hilbert space formulation of the relativistic quantum mechanical two-body problem was stud- 
ied by Rizov et al [5] . These authors derived an algorithm for a perturbative calculation of the 
energy eigenvalues in the case of a general energy dependent quasipotential. The problem was 
first solved for the particular case of the relativistic harmonic oscillator [6] . 

In recent years, numerous works have been devoted to the Hamiltonian formulation of 
relativistic quantum mechanics in connection with the manifestly covariant formalism with 
constraints [7, 8]. The main advantage of this approach is to allow the proper separation of 
relative and center-of-mass coordinates in few-body and even in many-body problem, and thus 
to eliminate spurious degrees of freedom. The constraints impose restrictions to the potentials 
to be physically acceptable. As a result, the potentials become energy dependent. Links with 
field theory have been developed by taking the Bethe-Salpctcr equation as a starting point. 
This makes the approach a challenging one. A typical application of the manifestly covariant 
formalism is a calculation of spectra of baryonic atoms. Work along this line has been carried 
out by Mourad and Sazdjian [9]. In particular, they have derived two-fermion relativistic wave 
equations, which acquire the form of Pauli-Schrodinger equations. 

The presence of the energy dependent potential in a wave equation has several non-trivial 
implications. The most obvious one is the modification of the scalar product, necessary to 
ensure the conservation of the norm [10]. However, the modification of the scalar product 
itself is not sufficient to justify the use of the common rules of quantum mechanics. In the 
next section, we illustrate the aforementioned difficulties on a few simple models, admitting 
analytical solution. 

Of course, testing formulas on a particular model does not guarantee their general validity 
in more realistic, physical systems. In section 3, we therefore analyze formal aspects and 
demonstrate under what circumstances a wave equation with an energy dependent potential 
corresponds to an acceptable quantum theory. We show that in such a case the wave equation 
can be transformed into an equivalent ordinary Schrodinger equation with a local, momentum 
dependent potential. Implications of the formal treatment are illustrated in section 4 and 
conclusions are drawn in the last section. 
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2 A toy model. 



2.1 Generalities. 

For the sake of simplicity, the toy model is constructed in the 1-dimensional space. ^ 
The starting point is the time dependent wave equation^ 

where V{x,i-^) denotes a real function of 2 variables. Setting ^{x,t) — e~*-^V(^) yields 

The first modification of the usual rules of quantum mechanics concerns the scalar product. 
This question has been studied by many authors in the one- and two-particle relativistic quan- 
tum mechanics (sec Sazdjian [7] and references therein). For completeness, we present a brief 
sketch in the 1-dimensional case : 

Consider two solutions of energies E and E' 

= e-'(^-'^)*$(a;), and ^,{x,t) = e-'^^^'-'^^'^ix), (3) 
with e — > 0. The usual continuity equation reads 

-P=--1 (4) 
dt dx^ ' ^ ' 

where 

p^<b:{x,t)^,{x,t) ; i = 

In the case of energy dependent potential, in order to fulfill the continuity equation, a term 

4n - nV:ix,t)[V{x,E) - V{x,E')]<^,{x,t) (5) 

at 



■^Throughout this section, we assume that solving the wave equation with an energy dependent potential 
we are able to select the states in such a way that each eigenvalue corresponds to a single normalizable wave 
function. This is not guarantee a priori, as will be discussed in the next section. 

^In this work we use m = % = 1. 
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has to be added to the left hand side of Eq. (4). After integration, Eq. (5) yields 

as e — > 0. As a result, in the hmit E' — > E, the scalar product (the norm) is given by 

N = J^^^*{x)[l-^Mx)dx. (7) 

Let us specify the stationary states by a global quantum number n. The orthogonality relation 
between two states n and n', n ^ n', is given by 

*n'W[l ^ ^ I'^nixjdx^O . (8) 

The generalization to higher space dimensions is straightforward. Expectation values and tran- 
sition matrix elements have to be calculated in a similar way, i.e. by including the correction 
term. Use has been made here of the continuity equation. The form of the scalar product 
with energy dependent potentials was also derived from Green's functions of quantum field 
theory [11]. 

The modification of the scalar product is a necessary but not sufficient step. A coherent 
theory requires the norm and the average values of positive definite operators to be positive. 
This is not guaranteed a priori. 

The second difficulty to face is the loss of the usual completeness relation X^n ^n(2^')^*t('^) — 
5{x—x'). Of course, this is a direct consequence of the fact that the functions \E'n(a:) constructed 
by the described procedure do not represent eigenfunctions of the same (linear selfadjoint) 
operator on L^(— oo, oo). It is tempting to propose : 

^*„(a;')*:(a;)[l - (^„„(x)] = 8{x - x') , (9) 

n 

V{x,Er,)-V{x,En) „ , , , , . dV{x,En) 
where ipnn'[x) = for n / n , and ipnn[x) = — . 

According to this ansatz, we have 

J^j{x)6{x-x')dx = <ilj{x') (10) 

= *„(x') / ■^j{x)-^l{x)[l-ipnn{x)]dx 

= -^jix') - J2 "^nix') I ^ j{x)^l{x)[ipnn{.x) - ifnAx)]dx . 
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The second term is non zero in general. Thus, the ansatz (9) is not vahd except for the case of 
a hnear energy dependence, when the terms (p{x)nj = ^{x) are independent of the state. 

The third comphcation arises due to the presence of the eigenvalue in the Hamiltonian 
(2). When calculating a commutator, one has to replace the eigenvalue by the corresponding 
operator. For the sake of illustration, let us take the following Hamiltonian : 



H^^ + Vo{x) + V,{x)E . 



(11) 



The commutator 



H,x 



is given by 



H, X 



-i(l-V,{x)) 'p. 



—ip + Vi{x) H, X 

It is to be stress here that P — —i-§^ not necessarily represents the momentum operator. 
Similarly, the double commutator takes the form 



(12) 



H,x 



,x 



^-{i-v,ix)y\ 



(13) 



In general, the commutators are more difficult to calculate, but the procedure is the same. It 
starts with 

• ^ n . dV{H,x)dH 



2.2 Explicit examples. 



In order to construct a toy model admitting analytical solutions, we consider the harmonic 
oscillator potential with an energy dependent frequency. In 1 dimension, the corresponding 
wave equation reads {cu — 1) 



(15) 



where 7 is a parameter (not necessarily small though most of our investigations are dedicated 
to small 7). Setting 

Xl^l + ^f{E^) , kl^2E^, (16) 



leads to 



dx"^ 



;i7) 
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The solutions are obtained by the ansatz 

= Cr,K{^|Kx)e-^-'''l^ , (18) 

which yields for \/A^^ = ^ 

^K{z) - 2z^K{z) + (^ _ i)/,^(^) = . (19) 
az^ az An 

Apart from the fact that A„ depends on the state, the solutions of this equation are formally 
the same as the Hermite polynomials. Their orthogonality is ensured by the weight function 

e-(^n+V)^V2[i _ 1 /(KQ - f{En) ^2] _ (^20) 

2 {En' — En) 

The energies and the A„'s are obtained from Eq. (16) together with 

kl^Xn{2n + l)^2En . (21) 
Two kinds of energy dependence have been considered : the linear and the square root cases. 



2.2.1 Linecir E-dependence 



For f{En) = En, — = and the conditions to be satisfied lead to a quadratic equation. 
Its general solution is given by 

The solutions have two branches, one of positive and one of negative values. However, only the 
positive energies are retained, since the negative ones are not normalizable. 

In Fig. 1, the energy En is presented as function of n for several values of 7 and n < 8. 
While positive values of 7 produce an increase of the level spacing as function of n, negative 
values of 7 compress the spectrum. For 7 < 0, is a monotonically increasing function of n. 
Asymptotically {n ^ 00), its derivative reaches — (2,^+1)3^3 ; ^-^d En ~ — ^. This limit ensures 
that A„ is real Vn, with A„ 0. 

As stated in the preceding section, physically acceptable potentials must yield positive 
norms. In the present case, we have 

Cl r e-^-^\l{^nx)[l - I x^]dx = 1 , (23) 

J —00 z 
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which gives 



^2 ^ 1 VKl 



4A 



1 1 -1 

' (2n + l) 



(24) 



Prom Eqs. (16) and (21), we get 

(2n + 1) 



7 = -^r; 



'1 + 



(2n + 1)2 



-r 



4 ' "V 16 
Consequently, the norm is positive definite Vn, irrespective of the sign of 7. 



(25) 



The modification of the scalar product affects also the calculation of expectation values. In 
the case of the mean square (ms) radius we obtain 



(2n + 1) 
2A„ 



1 - 



1_ 
4A. 



(2n + l) 



1 - 



37 [(2n + If + 1] 



8A„ (2n + 1) 



(26) 



For 7 positive, the sign of {x'^)n is finked to the sign of the last term. It is easy to show that for 
n — > 00 it takes the value 1/4, whereas for n = 0, it becomes negative beyond 7 = ^16/3. Thus 
the question is, for which state the ms radius is positive definite. For n = 1, it is straightforward 
to show (using Eqs. (16), (21)) that the last factor of Eq. (26) has no real root, and the ms 
radius is thus positive definite. 



As the critical value 7 = y 16/3 is large compared to 1, one might expect to avoid contra- 
dictions for 7 << 1. Unfortunately, explicit calculations of higher moments {x'')n {k even) for 
few low n states indicate that (up to the leading order): 



{x'')n OC [1 



7 {k 



4A. 



(27) 



Thus, for fixed n and positive 7, there will always be a value of k, beyond which {x^)n becomes 
negative. 

The above example clearly demonstrates that ad-hoc modification of the operation which 
should play the role of the scalar product and the positivity of the "norm" of the assumed 
energy eigenfunctions are not sufficient to ensure the positivity of the expectation value of a 
positive definite operator. It underlines an internal inconsistency for 7 > 0. 

As far as the closure relation is concerned, as stated in the preceding section, the ansatz (9) 
seems to be valid for a linear energy dependence. It is straightforward to get a numerical test 
by calculating, for instance, 

{x\ = j:\{n\x\0)\' . (28) 
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Up to linear terms in 7, we verified that 



{l\x 



(29) 



Higher contributions, |(n|a;|0)p , where n > 3, are proportional at least to 7^. This results 
from the fact that in the ordinary harmonic oscillator (7 = 0), the x operator connects only 
states which differ by one quantum of energy. It is obviously a tedious work to check Eq. 
(28) to all orders in 7. However, since we know the exact value of (x^)o, we can estimate the 
contributions from the first excited states and discuss the convergence of the series expansion. 
As an example, we have calculated the first two terms for 0.01 < I7I < 0.50. The results are 
displayed in Table 1. They are compared with the exact values and with the estimates obtained 
keeping only contributions linear in 7. Whereas the latter approximation becomes inaccurate 
by few percent beyond 7 = 0.25, the convergence seems rapid. The contribution of the third 
excited state is less than 2 % for 7 = 0.5. 



As the second example, we have tested the same way 



+ 2b^ 



4(1 - , 



(30) 



keeping only terms linear in 7. In Table 2, we compare the exact values with the linear estimates 
and with the sums of the series expansion as a function of Umax (the highest state involved in 
the sum). The situation is somewhat worse than for < >o : the hnear approximation is less 
efficient and the convergence of the series expansion is slower. 

The above results support the extension (9) of the closure relation. They also underline the 
care one has to take when using 7 as a small parameter in a series expansion. On the other 
hand, for 7 > 0, we note a tendency to overpass the exact value, which appears beyond 7 = 0.4 
for < >o and beyond 7 = 0.35 for < >q. We have verified that this effect is not due to 
numerical uncertainties. Again, it is a signature of the deficiency of the algorithm for 7 > 0. 

The consequence of the modification of commutation relations can be illustrated on the 
dipole sum rule. We recall that it is obtained by averaging over the ground state the double 
commutator [H, x\,x . 



According to Eq. (13), we have 



[H, x\,x 



and the ground state expectation value yields 



^(K-^o)|(n|a;|0)p 



1 



1-7x72 



[H.xlx |0) 



-] 



[1- ' 



(31) 



(32) 
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As a numerical test, we adopted several maximum values rimax of n considered in the sum. 
Results are displayed in Table 3 for rimax = 1 up to 7. They are compared with the exact value 
given by the last term of Eq. (32). The convergence rate is slower as 7 increases, as expected. 
The convergence towards the right result for the low 7 values supports the validity of the rules 
used to derive Eq. (32). For 7 > 0, contrary to the < x'^ >o and < x"^ >o cases, the sum seems 
to converge toward the exact value. This comes from the fact that the leading contribution to 
the ground state average of Eq. (32) is the norm, for which the orthonormality of the wave 
functions ensures that J2n I < n\0 > |^ = 1- 



2.2.2 ^/E - dependence 



The case is chosen for the sake of comparison with the results of the previous example. 

A smoother energy dependence should generate weaker effects. On the other hand, the simpler 
generalization (9) of the closure relation is valid only approximately. Therefore, we do not 
expect to calculate expectation values of commutators (13,14) with unlimited accuracy. 

The eigenvalues of Eq. (15) are in this case given by 

4El-{2n+lf[l + j^n]^0 . (33) 

On the real half-plane > 0, for each sign of 7 and each n, there exists a single real eigenvalue 
leading to square integrable wave function. The spectrum exhibits the same features as the 
linear E'^-dependence. 

For small positive 7 (typically 7 < 0.2) and few lowest states (n < 3), a perturbative 
estimate yields 

K-(n+^)[l + |^^^^]. (34) 



For 7 < 0, the large n limit is given by 



hm E^^^. (35) 

n— >0O ry-i 



As in the linear £^„-dependence, it ensures A„ to be real and reaching as n — >^ 00. 



The norms and the (assumed) predictions of the average values of observables are calculated 
by introducing the multiplicative factor 

1 - ——=^^—=-x'^ . (36) 

2(V-E'n + vEn') 
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The norms, given by 



Cl 



1 



1 (2n + 1) 



'Er, 



(37) 



are positive definite as can be proved using Eq. (21) and (33). 
The ms radius 



2n + l 

2Xn 



1 - 



1 7(2n + l) 



1 - 



3 7[(2n+ 1)2 + 1] 



16/E;: A„(2n + 1) 



(38) 



is positive definite Vn, as well. On the other hand, for {x^)ni the same reasoning as for the 
linear case yields (up to the leading order) 



{X^)n OC [1 - 



4A„v^ 



(39) 



For fixed 7 > and n, there exists a value of k beyond which {x^)n becomes negative. In spite 
of a smoother energy dependence than in the linear case, very similar problems occur in the 
case of expectation values of positive definite operators for 7 > 0. 

Because are energy dependent, the "closure ansatz" has no chance to be given by the 
simple extension (9). Nevertheless, one can estimate the correction terms. To fix the order of 
magnitude, we take Eq. (10) for = ^0, where 



7 X 



2\/En 



'Eq + V%i) 



(40) 



We consider the first two contributions to the sum of Eq (10). Because of parity, they arise 
from '^2 and \E'4. The results are displayed in Table 4. The correction is small in the range of 
considered 7 values, but not necessary negligible. It is expected to decrease with increasing n. 

Neglecting the correction terms, we can perform the same tests as in the preceding subsec- 
tion. Up to linear terms in 7, 

1 7 1 
(a;2)o~-(l 1=) is to be compared to |(l|a;|0)|2 a; -(1 - O.657) . 

2 -\/2 2 

The difference between the two estimates is of the same order as the correction terms listed in 
Table 4. 



A more quantitative estimate is displayed in Table 5. The exact value of {x'^)q is compared 
with the sums Xln""""" with rimax = 1,2 and 3. The series expansion converges rather 

rapidly but not to the exact value. The discrepancies are of the same order as the one observed 
in the above expression limited to terms linear in 7. 
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2.2.3 Further energy dependences 



The above two examples yield an overview of the difficulties linked to the energy dependence 
of the potential in a wave equation. There is no need to pursue this investigation much further. 
One point, however, is worth mentioning. It concerns the cases when higher powers of the 
energy dependence are considered. Take for instance /(£^„) = E^. It is easy to verify that the 
positive eigenvalues are given by 

E.^^±^. (41) 

We note that for positive 7 there will always be a critical n beyond which the eigenvalue 
becomes complex. 



3 A quantum theory. 



A toy model is well suited for illustrative purposes. However, in order to go beyond and check if 
it corresponds to a well defined theory, a more fundamental analysis is required. Mathematical 
aspects of theories with energy dependent potentials have been studied in ref. [12,13,5]. In 
particular, these works are devoted to the positivity of the norm and the occurance of patholo- 
gies [1]. 

Let us define (for a given Hilbert space H) a Hamiltonian 

H{z) =Ho + V{z) , (42) 

where Hq is a selfadjoint operator and there exists such a set Ai dTZ that for any fixed value 
z & A4 each of the operators V{z), H{z) is selfadjoint too. 

Let us make attempt to construct a quantum theory (QT) in such a way that 

1) the vectors |^ >G Ti., at least some of them, correspond to states of a physical system 
(here and in the following "state" is used in the sense of "pure state" only), 

2) the time evolution of these states is determined by the equation 

^|l*W>-^(^|)l*W> ■ (43) 

First of all, one should keep in mind the following requirements, which have to be satisfied by 
any QT : 
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i) Each state of the considered physical system is mapped to a ray in the corresponding 
Hilbert space (Tip). 

ii) Each state of the considered system represents a hnear combination of its stationary states. 

iii) The initial state of the system (its state at to) uniquely determines its states at each later 
moment t > (provided the system was not submitted to any measurement and/or observation 
during the interval {to,t > ). 

The last requirement apparently limits the acceptable 2;-dependence of V{z) to the linear 
one only. But one can try to avoid this sharp restriction in the following way : the Hilbert 
space Hp need not be identical with H. In fact, the two spaces might differ in 2 aspects : 

a) some of the vectors from H need not belong to Tip] these elements of 7i would not describe 
any state of the considered system, i.e., the vector space Tip could be a (proper) subspace of 
the vector space H. 

b) The vector space Tip could be endowed with a different scalar product (( ■ | ■ > ) than 
the vector space Ti. (We reserve the symbol < ■ | ■ > for the scalar product in the Hilbert space 
H). In other words, the Hilbert space Hp need not be a subspace of the Hilbert space H. 

From the requirement i) directly follows that a vector 

I* >e Hp , I* >7^ (44) 

corresponds to a stationary state of the considered system if and only if there is such a real 
number E that the vector valued function (of time) 

\qf(t) >= 6-^-^*1* > (45) 

satisfies Eq. (43). Moreover, it is evident that the function (45) solves Eq. (43) if and only if 
the vector (44) satisfies the relation 

(^0 + V{E) ) I* >= £; I* > . (46) 

It is customary to call this relation the (time independent) "Schrodinger equation with an 
energy dependent potential". 

Let us recall that \/E e M there 3 >eH; 

{Ho + V{E) ) I* > e H . (47) 
On the other hand (usually) for most of values E & M. none of 

I* > e 7i , I* >v^0 (48) 
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satisfies the condition (46). To find out those values of E for which Eq. (46) has nontrivial 
solution in 7i, let us first of all consider the equation 

H{z)\'^ >= E{z) I* > . (49) 

For simplicity, we will limit our discussion to the case where & A4 the spectrum of H{z) is 
purely discrete and non-degenerate. This means that for any fixed value z e Al, to each 
eigenvalue E^iz) of H{z), there 3 \'${z) > eH : 

H{z)\^n{z) >= E^{Z) \^n{z) > , \^n{z) >^ (50) 

and no other linearly independent solution of Eq. (46) with E = En{z) exists in Ti. It may 
happen that for a given value of n there exist such values z & A4 that the relation 

En{z) = z (51) 

holds. In the following, we will enumerate by the index n only those eigenvalues of H{z) for 
which Eq. (51) has at least one solution in JH. Those values of z E A4 which satisfy the 
relation (51) will be denoted as 

, m=l,...,M„>l, (52) 
i.e., z^^ are those elements of M. for which 

Eniz^r^^) = Zt^ . (53) 

Due to Eq. (50) we have 

Hzt^) \'^n{zt^) >= E^{zt^) 1^.(4")) >= 4') |v]/„(4")) > . (54) 

Moreover, due to the assumed non-degeneracy of the H{z) spectrum, we know that for each 
given z'fil^ the last equation has one and only one linearly independent solution in Ti. Let us 
denote it as |^'^-' >. Without any loss of generality, on can impose the normalization condition 

< I >= 1 . (55) 

In order to simplify symbolics, let us switch from the two-index to one-index notation 

n-l 

{(n), m} — J = Mk + m , where Mq = (56) 

fe=0 

and 

^,■ = 4"), \^^>= |*W> . (57) 
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If one accepts the "philosophy of energy dependent potential" , then each of the values 

, j^l,...,N = ^Mk (58) 

k 

represents an energy level of the considered system (and each energy level coincides with one 
of the values (58)) and the vector > unambiguously describes the corresponding stationary 
state. 

In order to satisfy the aforementioned condition ii), one has to require Vj^' >e Hp to be 
expressible in the form 

1^ >= ^ Cj l^j > , where cj — const. . (59) 



Such a vector space Tip could represent Hilbert space assigned to a physical system in the 
framework of any QT only if it could be endowed with such a scalar product ( ■ | ■ > that the 
orthogonality condition 

(*,|^'fe>=o , yj^k (60) 

would be guaranteed. On the other hand, 

{Ho + V{z,) I ^j>=z, |*,> (61) 
implies (due to selfadjointness of Hq and V{z)) that the identity 

{zj - Zk) < >=< -^kliViz,) - V{zk)]\-^, > (62) 

is valid Wj k. 

Hence, the requirement (60) would be automatically satisfied if the scalar product in Hp is 
defined in such a way that 

(*,|*fc>=<*,|iy,fe|*ifc> , (63) 

where 

Wj, = 1 - Qjk (64) 

and 

4,^ ^fe)-^(^^) , Vj#^ . (65) 

Zj Z]f 

In order to define the scalar product in Hp, it remains to specify the meaning of the r.h.s. 
of the expression (65) for j — k. It is tempting to take 
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We are not going to dwell on arguments neither in favor nor against this specific choice of the 
"scalar product" ( ■ |- >, utilized e.g. in Eq. (8). Let us stress instead that one can face serious 
problems when introducing any scalar product ( • |- >, which would guarantee fulfillment of 
the requirement (60). The point is that the conditions (54) and (55) define the set of vectors 

|*, > , j = l,...,N (67) 

(up to phase factors) uniquely but they do not imply the linear independence of all the vectors 
(68). One should not forget that the requirements (54) and (55) do not guarantee that all the 
vectors (67) are eigenvectors of a common selfadjoint operator. On the other hand, if 

then from the condition (60) follows 

(v^,|v^,- >= , (69) 

i.e., l^^j > would be the null vector. But from the condition (55) we know that it is not so. This 
means that if not all the vectors (68) are linearly independent, than it is impossible to define 
any scalar product ( ■ |- > in Hp in such a way that the condition (60) is satisfied. Especially, 
this means that in such a case 

i) the relations (63)-(66) do not define a scalar product 

ii) it is impossible to modify the r.h.s. of equation (66) in such a way that the relations 
(63)- (66) would define a scalar product in Tip. 

In other words, in such a case, the aforementioned algorithm does not represent any QT. 

On the other hand, it is clear that for any given quantities ^o, there exists (at least 
one) such a set Ai that all the vectors (67) are linearly independent. Therefore, at least in 
principle, one can avoid the above indicated problem simply by a "proper choice" of the set 
M. 

Of course one should keep in mind that 
a) the requirement of the linear independence of the vectors (67) does not determine the set 
A4 uniquely. ( It could be tempting to choose Ai as the largest subset of real values z for which 
the operators H{z), V{z) are selfadjoint and all the vectors (67) are linearly independent. But 
generally there still could be more then one subset of TZ which fulfills all these conditions.) 
/?) the "energy spectrum" produced by such an algorithm substantially depends on the par- 
ticular choice of the set M.. 

The main moral from these observations is that the specification of the operator (valued 
function) (42) itself need not determine uniquely the energy spectrum of the corresponding 
physical system with energy dependent potential. 
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To proceed further in our analysis, let us from now on assume that the set Ai is already 
chosen in such a way that all the vectors (67) are linearly independent, and therefore the vector 
space Hp (59) is a A^-dimensional subspace of the vector space H. Let us also define vectors 

l^i ^0) 1^^. > , (70) 

where A^^^ are non-vanishing constants determined (up to their phase factors) by the require- 
ment 

(vl>i|vl/^'>=<5^., ; j,A; = l,...,7V, (71) 
once the scalar product ( • | • > in Tip is specified. 

Since Hp is an A'"-dimensional subspace of H, there exist such vectors 

\Xj>eHp, j = l,...,7V (72) 

that 

< Xj\Xk>^ Sjk , j,k ^1,...,N , (73) 

and any vector |^ > e Hp can be expressed as a linear combination of them. Especially each 
vector (70) is expressible as 

N 

\^'>=J2^kj\Xk>,J = l,...,N (74) 

k=l 

and the constants 

Cfe,- =< Xfel*'' > ■ (75) 

Evidently these relations could be expressed also in the form 

|*^ >=Clx, >, J = l,...,N , (76) 

where 

TV N 

i = Y: I*' X Xk\ = E ^eklxe X Xk\ (77) 

k=l l,k=l 

is a linear operator (on Hp) and 

Ckj =< XklilXj > ■ (78) 

Moreover, due to the linear independence of all N vectors (70), the relations (76) have to 
be invertible, i.e., the operator ^ is nonsingular. Therefore the relations (76) are equivalently 
expressed as 

IXi >= ?7l*'' > , j^l,...,N , (79) 
where the linear operator (on Hp) 

V = ■ (80) 
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Let us define V |^ > e Tip the corresponding vector (in Tip) 

\^ >= > . (81) 

Recalling Eq. (79) one can see that 

>- \Xj > (82) 

and hence the equivalence 

N N 

1^ >= $:c,|^^- >4^\^>= J2cj\Xj > (83) 

holds. Therefore, for any vector >, |$ > G Tip the relation 

($1^ >=< $||r > (84) 

holds. It means that the scalar product ( • | • > could be expressed in terms of the scalar product 
< - j- > as 

($1^ > = < $|«;|^ > , (85) 

where the linear operator 

k = rj'^rj . (86) 

It is useful to express Eq. (85) as a relation between two different bra-vectors related to the 
same ket by the two different definitions of the scalar product (on the same vector space Hp) : 

($1 =< . (87) 



One should not overlook that using this formalism one can e.g. express the closure relation 
corresponding to the orthonormality condition (71) in the following simple form 

N 

^I^J' > (^^1 = 1 . 87a) 
Needless to say that this relation is equivalent to the closure 

N 

J2\Xj ><Xj\ = 1 ■ 



In order that the considered algorithm could provide a QT of a physical system, one has 
to assign to each observable A the corresponding sclfadjoint operator A. At this point, it is 
important to keep in mind that the notion of conjugation of linear operators depends on 
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the choice of the scalar product. We have tacitly used the symbol to denote the operator 
conjugated to the operator A in the sense of the scalar product < - I- >, i.e., the relation 

< >=< >* (88) 

holds. 

Let us introduce the symbol A* to denote the operator conjugated to A in the sense of the 
scalar product ( • |- >, i.e.. the relation 

(*|i#|$ >= ($|i|* >* (89) 

holds. By using the relation (87) one finds 

($|i|* >*=< ^\kA\^ >*=< ^\A+k\^ >= {^\k-^ A+ k\^ > (90) 

and therefore 

A* = k-^ A+ k = 11+ i+ fj+r] . (91) 
One can utilize the formal relations 

(I* > )+ ^1 ; (I* > )# = (*| , (92) 

and therefore 

(l^r X 01)+ = 10 X *| ; (I* > (0|)# = 10 > (*| , (93) 

etc. It is also useful to notice that the operator (86) is selfadjoint from the point of view of 
both scalar products, i.e. 

k^k+ ^k* . (94) 

Needless to say that when we required (below Eq. (87)) the operator A, assigned in the 
considered algorithm to an observable ^4, to be selfadjoint, we meant that the relation 

A*^A (95) 

has to hold. 

The discussed algorithm may appear as a sort of QT which is still "pecuhar" in two aspects: 

i) the (total) energy is the only observable to which no selfadjoint operator was assigned, 

ii) evolution of the state has not been described by an "ordinary Schrodinger equation" . 

Let us show now that both of these points could be circumvented, namely, that the consid- 
ered algorithm may be equivalently expressed in the form of an ordinary quantum 
mechanics: 
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Let us define the operator 

TV 

^ = ^ZjI*'' > (*^| . (96) 

From Eqs. (93), (54) and (70) follows that the relations 

H*^H , (97) 

^l^i >= Ejl"^^ > , where Ej = Zj , (98) 

are guaranteed. Moreover, from Eqs. (45), (59) and (70) follows that the solution (43) with 
the initial condition 

l^(to) >= I* > e (99) 



may be expressed in the form 



-^{t) >= ^dj-e"^^^*!*^' > , where dj = (100) 



and hence the "ordinary Schrodinger equation" 
holds. 



i|-|^(i) >=^|^(i) > (101) 



Therefore the considered QT (with an "energy dependent potential" ) could be equivalently 
rephrased as an "ordinary quantum mechanics" (QMl) in such a way that 

1) the Hilbert space of the corresponding system is identified with the vector space Tip C H. 
endowed with the scalar product ( • |- >, 

2) the Hamiltonian is identified with the operator (96). 

Actually, one can rephrase the considered QT also as an "ordinary quantum mechanics" 
(QM2) of the same physical system in such a way that the Hilbert space assigned to this system 
is identified with the subspace Tip C 7Y, i.e. with the vector space Tip endowed with the scalar 
product < - j- >. 

To this end, it is sufficient to realize that once one associates to each linear operator A the 
operator 

A = i^Ai , (102) 
then, from the definition (81), one immediately sees that V|^ > G Tip 
the relation |$ >= A > is equivalent to |<l >= A\^ >, 
and the relation A* = B might be equivalently expressed as ^4 — B. 
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On the other hand, from the formulae (96), (87) and (82), one immediately obtains 



N N 

H^Y. ^i^l*' >< = E ^j\Xj >< Xj\ ■ (103) 

j=l j=l 

In this case, the stationary states are described by the vectors \xj > and the evolution of 
the states of the considered system is described by the Schrodinger equation 

ii^(i) >= H\^{t) > . (104) 
The operators associated to observables are now selfadjoint in the sense of the relation 

1^ = 1 . (105) 

We have seen that the considered algorithm could represent a QT if and only if the operation 
( ■ |- > (introduced as a part of this algorithm) is expressible in the form (85). where «; is a 
selfadjoint positive definite operator. How can one check whether this condition is met? 
First, it would be useful to express the operator fj in the form 

TV 

^ = E ^jklXj >< Xk\ ■ (106) 
j,k=i 

One finds that 

Vjk =< xMXk >=< *^|r/+r/|xfc >= (^^Ixfc > • (107) 

This means that the coefficients on the r.h.s of Eq. (106) are identical with the coefficients in 
the expression 

N 

\Xj>=Evk,\'^'> ■ (108) 

k=l 

Therefore one can obtain the operator corresponding to any (arbitrarily chosen) orthonormal 
basis (in the sense of relation (73)) 

{ \xj>, j = l,...,N} (109) 

where each of its elements is expressed as a linear combination (108) with the coefficients rjkj 
utilized in Eq. (106). 

Once the operator rj is known, one can check whether all the eigenvalues of the operator 
(94) 

k = fj+fj (110) 
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are non- vanishing. (Needless to say that the operator k is uniquely determined by the operation 
( ■ |- >. The freedom in the choice of the operator fj is just a trivial consequence of the 
fact that predictions of any (ordinary) quantum theory are invariant under arbitrary unitary 
transformations.) If it is the case then the considered algorithm could represent a QT. In 
the following we shall already tacitly assume that the operation ( ■ |- > is introduced in the 
considered algorithm in such a way that it can be interpreted as a scalar product in Tip. 

We have seen that this QT might be equivalently rephrased as a QMl as well as a QM2. 

For passing from the QT to the QM2, the knowledge of the operator fj appears essential (see 
Eqs. (102), (103) and (81)). We have already discussed the steps which could in principle lead 
to the construction of fj. Unfortunately, to proceed this way in practice usually represents a 
formidably complicated task. One naturally asks whether there exist any simpler way leading 
to the same goal. 



In fact, it is easy to find a more straightforward procedure, at least for one particular case. 

Let 

V{z) = zV (111) 

and the operator ( • |- > is defined by the formulae (63)-(66). In the case (111), the formula 
(64) simplifies to 

Wjk^l-V , (112) 

and therefore (cf. Eqs. (63) and (85)) 

k=l-V . (113) 

Such an algorithm might represent a QT only if all values from the spectrum of the operator 
V are smaller than 1. The evolution equation (43) can be rewritten as 

ij^il-V)\^{t)>^Ho\^it)> (114) 

or equivalently as 

if^\^{t)>^iHom{t)> : (115) 
where we have used the notation (81): 

|^(t) >=7/|^(t) > , (116) 
and without any loss of generality we have identified 

fl^f}+ = Vk^\/l-V. (117) 

The equation (115) is of course nothing else than the (ordinary) Schrodinger eq. (104) with 
the Hamiltonian 

H = ^^^0^^ , (118) 
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where 

i=[l- V]-'/' . (119) 
Therefore, in this case we have reached the QM2, which is equivalent to the QT, just for free. 

Let us recall that in the framework of QM2, to each observable A is assigned an operator A, 
which fulfills the condition (105). In particular (in the considered case), the operator assigned 
to the (total) energy is given by formula (118). 

From Eq.(102) follows that the operator A, which is in the framework of QMl associated 

to the observable A, is related to the operator A associated to the same observable in the 
framework of QM2 as follows 

A = ^Afj , 

i.e. 



[l-V) 

Especially, the Hamiltonian in the QMl is the operator 



A^/{l-V). (120) 



I.e. 



H^CHfj^eHo, (121) 



H = ^Ho . (122) 

l-V 



4 The toy model revisited. 

In the light of the preceding section, we are going to re-examine our toy model. We will present 
an example of QT which reduces to QMl or QM2. Let us consider the special case when 

n = L^{-oo,oo) , (123) 

^o^-^^ + W, (124) 

V = V{x) , (125) 
and the real function V fulfills the condition 

V{x) < 1 Vx e (—00,00) . 

Physically, this should be a one dimensional version of QT of a particle in an external field. 
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In this case 



To be more specific, let us take 



rf2 V d I 2V"{l-V) + W 
dx^^ 1-V dx^4 (1 - 1/)2 



V{x) = -Kx^ , K>0 
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Then 



1 



2Kx d 



— -K 



1 - 2Kx'^ 



(126) 



(127) 



(128) 



-dx'^ 1 + Kx^ dx (1 + Kx2)2- 

To obtain the energy spectrum, one has to find those values of the parameter E for which the 
differential equation 



{-I 



d^ 



2Kx d 



— -K 



1 - 2Kx'^ 



+ V^{x)}xe{x) = E{1 + Kx^)xe{x) (129) 



Lda;2 l + Kx"^ dx {I + Kx'^f 
has solutions in Tip of QM2, i.e., which satisfy the integrability condition 

/oo 
\XE{x)\^dx < oo . 
-oo 



(130) 



(For simplicity sake, we are explicitly describing here the procedure which concerns only the 
discrete part of the energy spectrum). 



To be even more specific, let us consider the special case when 

Vo{x) ^Ax^, A > . 

Then Eq. (129) yields 



1 „ Kx 

2^ ^ (1 + Kx'') 



X- 



K 1 - 2Kx'^ 



2 {l + Kx'^y 

In the asymptotic region oo, Eq. (137) acquires the form 

where 

\^{E) = 2{A~KE) . 



+ E + {EK - A)x^]x = 



(131) 



(132) 



(133) 



(134) 



Thus, for each real value of the parameter E, one can choose 2 independent solutions of Eq. 
(132) in such a way that their behavior in the asymptotic region a; — oo is dominated by the 
factor 

(135) 



HE) ^2 



HE) ^2 



e 2 ^ and e 2 , 
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respectively. Similarly, one can choose 2 independent solutions of (132) in such a way that their 
behavior in the asymptotic region x — > —oo is also dominated by the factor (135). 



For E > ^, X{E) is purely imaginary, and consequently both independent solutions oscillate 
in both asymptotic regions. This is why we strongly suspect that each value E E (^,oo) 
represents a (twice degenerated) value from the continuous part of the spectrum of the operator 

H. We are not going to dwell on closer examination of this point here. 

Let us turn to the examination of Eq. (132) for < ^. Due to the invariance of the 

operator H under the replacement x — > —x, for each value of the parameter E < A/K, there 
could be at most one solution of Eq. (132) belonging to L^(— oo,oo). This means that the 
whole discrete part of the spectrum is non-degenerated. To find the explicit formula for this 
part of the energy spectrum, let us express the function Xe{x) as 



X{x) = ^JX + Ky^fiy)e-y , where y = VXx . (136) 
Then Eq. (132) implies that the function / has to solve the equation 

2E 

r-2?//'+(^-l)/ = 0. (137) 

It is well known from elementary QM that 
i) if 

2 E 

— = 2n + 1 , (138) 
A 

where n is an even (odd) nonnegative integer, then in the asymptotic region \y\ oo, symmetric 
(antisymmetric) solutions of Eq. (137) 

f{y) cx HM , (139) 
where is the Hermite polynomial, whereas antisymmetric (symmetric) solutions behave as 

f{y) e^' . (140) 



ii) If the condition (138) is not satisfied, then the general solution of Eq. (137) has the asymp- 
totic behavior (140). 



Therefore, the discrete part of the energy spectrum is indeed non-degenerate and its energy 
levels are determined by Eqs. (138) and (134) as 



E„ 



2n+l 



^JSA + {2n + IfK^ - {2n + 1)K 



(141) 
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One should not overlook that E^, is a monotonically increasing function of n with the upper 
limit 

lim K = 4 ■ (142) 

n— >oo 

The stationary state with the energy (141) is described in the framework of QM2 by the 
wave function 

Xn{x) = C„Vl + Kxm^{^^x)e-^-^"^^ , (143) 
where C„ is a normalization constant and 



Xn = 2^A-KEn. (144) 

In the framework of QMl, the same stationary state is described by the wave function 

r{x) = ixnix) = Cr,H,,{Jx'nx)e-^-^'/' . (145) 



The same wave function describes this state also in the QT. 



The above analysis sheds some light on the toy models of section 2. In particular, the linear 
^'-dependence with 7 < coincides with the example presented in this chapter by setting 
A — 1/2 and K — In this special case (as one can see from formulae (92), (118)) the ansatz 
(9) looks like relation (87a) expressed in the "x-representation". It means that this ansatz 
would indeed be valid provided that 

a) the "physical Hilbert space" Tip is isomorphic to L^(— cxo, 00) 

b) there is no continuous part of the energy spectrum. 

As we have already stressed these conditions are apparently not met. Therefore, even in this 
special case the r.h.s. of the formula (9) is only an approximation of its l.h.s. 

Positive values of 7 are unsuitable. According to Eq (117), the operator rj reduces to zero at 
some place, becomes complex and is not invertible, no matter how small is 7. This is the origin 
of the questionable results concerning the sign of < x*^ >„ or the overestimate of < >o and 
< >o when applying the closure relation (Tables 1 and 2). The situation is similar for the 
A/E^-dependence as regards to the sign of 7. However, constructing the corresponding rj and k 
operators can be quite tedious. It may also happen that a judicious choice of the \x > vectors 
suggests a useful approximation. It lies outside the scope of the present work to investigate 
this point further. 

Let us now return to the more general case of the one dimensional QT determined by Eqs. 
(123)-(125): the Hamiltonian (126) can be expressed as 

H^T + H^i), (146) 
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where the first term 

T=^ (147) 
describes the kinetic energy of the considered particle. 

1 f_ 1 r2\/"(l - 1/) + 3y'2 . V 



H 



{I) 



-^{Vo-- r ' r \" + - VP ] I (148) 



4(1 -\/)2 

describes its interaction with the external field and 

P = -i-^ (149) 
ax 

is the operator assigned to its momentum in the framework of QM2. 

It should be noticed that 

i) rephrasing the QT (with energy dependent potential) as QM2 (ordinary quantum mechanics) 
revealed that the interaction of the considered particle is momentum dependent. 

ii) The expression for the Hamiltonian H describing the considered physical system in the 
framework of QMl is obtained from Eqs. (146) - (149) simply by omitting all tildes in these 
expressions. 

iii) In the framework of QMl the momentum is described by the operator (cf. (102)) 

p = iPij = p+i[p,ii], (150) 

dx 2 ^1 _ v{x) 

iv) If one insists to work in the framework of QT with an energy dependent potential, one has 
to identify again the momentum operator with the operator (151). 



5 Conclusions. 



The present work is devoted to wave equations with energy dependent potentials. Calcula- 
tions performed in the framework of several toy models clearly demonstrated that "plausible" 
modification of the "scalar product" does not guarantee intrinsic consistency of the resulting 
" quantum theory with an energy dependent potential" . 

In order to show what sort of caution one has to observe towards such straightforward 
implementation of energy dependent potentials in quantum theory we have resorted to a more 
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fundamental analysis of this problem. The main moral following from the analysis could be 
summarized as follows: 

i) Proper care has to be paid to the assignment of Hilbert space to the considered physical 
system. Especially one has to observe that 

a) vectors corresponding to stationary states with different energies have to be orthogonal ^ 

b) closure relation expressed in terms of stationary state wave functions properly reflects their 
normalization 

c) operators corresponding to observables are all indeed selfadjoint. 

ii) Once a " quantum mechanics with an energy dependent potential" does not suffer by intrinsic 
inconsistencies one can rephrase it into a form of ordinary quantum mechanics. The price one 
has to be ready to pay for such a transition to "ordinary theory" is a possible (additional) 
momentum dependence at that part of the Hamiltonian which should describe interaction. 
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Figure 1: En s& function of n for the linear E-dependent harmonic oscillator, with 7 = 
0, ±0.1, and ± 0.2. The solid hne corresponds to 7 = 0, points above (below) correspond 
to 7 > (7 < 0). 
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Table 1: Comparison of < >o as function of 7 for linear E dependence with predictions 
of the linear approximation (3rd column) and the sum (28); rimax denotes the highest state 
involved in the sum. See text for details. 



7 


< >o 


- §7) 




^max 3 


0.01 


0.496258 


0.496250 


0.496253 


0.496258 


0.05 


0.481445 


0.481250 


0.481338 


0.481445 


0.10 


0.463281 


0.462500 


0.462894 


0.463275 


0.25 


0.411124 


0.406250 


0.409449 


0.411012 


0.50 


0.331895 


0.312500 


0.329503 


0.332217 


- 0.01 


0.503758 


0.503750 


0.503753 


0.503758 


- 0.05 


0.518945 


0.518750 


0.518818 


0.518945 


- 0.10 


0.538281 


0.537500 


0.537731 


0.538270 


- 0.25 


0.598624 


0.593750 


0.594457 


0.598111 


- 0.50 


0.706895 


0.687500 


0.686108 


0.699289 
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Table 2: Comparison of < a;^ >o as function of 7 for linear E dependence with predictions of 
the linear approximation (3rd column) and the sum (30). See text for details. 



7 


< >o 


!(i-i7) 




^max 4 


^max 6 


0.01 


0.738797 


0.738750 


0.738760 


0.738797 


0.738797 


0.05 


0.694912 


0.693750 


0.694091 


0.694904 


0.694911 


0.10 


0.642105 


0.637500 


0.639288 


0.642006 


0.642097 


0.25 


0.496767 


0.468750 


0.487138 


0.495611 


0.496552 


0.50 


0.294496 


0.187500 


0.292966 


0.301638 


0.302978 


-0.01 


0.761297 


0.761250 


0.761259 


0.761297 


0.761297 


-0.05 


0.807432 


0.806250 


0.806378 


0.807421 


0.807432 


-0.10 


0.867270 


0.862500 


0.862603 


0.867066 


0.867248 


-0.25 


1.061827 


1.031250 


1.025211 


1.053279 


1.057984 


-0.50 


1.439879 


1.312500 


1.263341 


1.340323 


1.361776 
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Table 3: Test of the dipole sum rule for linear E dependence and 7 > for several maximum 
values Umax in the sum (32). The exact value is given by the last term in Eq. (32). 



7 


_ 

'^max ^ 




- ^ 

'^max "3 


_ 

^max ' 


exact 


0.01 


0.501236 


0.501250 


0.501250 


0.501250 


0.501250 


0.05 


0.505894 


0.506247 


0.506249 


0.506250 


0.506250 


0.10 


0.511061 


0.512451 


0.512492 


0.512495 


0.512496 


0.25 


0.522114 


0.529647 


0.530583 


0.530834 


0.531189 


0.50 


0.526582 


0.545979 


0.550165 


0.551709 


0.562017 


-0.01 


0.498736 


0.498750 


0.498750 


0.498750 


0.498750 


-0.05 


0.493404 


0.493748 


0.493750 


0.493750 


0.493750 


-0.10 


0.486139 


0.487461 


0.487500 


0.487503 


0.487504 


-0.25 


0.460802 


0.467450 


0.468276 


0.468497 


0.468811 


-0.50 


0.410368 


0.425484 


0.428746 


0.429950 


0.437983 



Table 4: Test of the "modified" closure relation (10). ^'2 and ^'4 denotes the 1st and the 2nd 
contribution to the sum (10), respectively. 





7 > 


7 < 


I7I 




^'4 


^2 


^'4 


0.01 


-0.000428 


0.000000 


0.000426 


0.000000 


0.05 


-0.002122 


-0.000030 


0.002145 


-0.000033 


0.10 


-0.004211 


-0.000116 


0.004301 


-0.000139 


0.25 


-0.010212 


-0.000628 


0.010670 


-0.000954 


0.50 


-0.019131 


-0.001965 


0.019540 


-0.003851 
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Table 5: < x"^ >o as function of 7, 7 > 0, for \/^ dependence. Exact values (exact) are 
compared with the sums X^^""""' | < n|a;|0 > rimax = 1,3 and 5. See text for details. 



7 








exact 


0.01 


0.498530 


0.498534 


0.498534 


0.498229 


0.05 


0.484265 


0.484270 


0.484280 


0.482782 


0.10 


0.469378 


0.469396 


0.469434 


0.466445 


0.25 


0.429326 


0.429411 


0.429591 


0.422212 


0.50 


0.375328 


0.375554 


0.376051 


0.361802 


-0.01 


0.501465 


0.501469 


0.501469 


0.501765 


-0.05 


0.516501 


0.516507 


0.516521 


0.518156 


-0.10 


0.533654 


0.533682 


0.533745 


0.537310 


-0.25 


0.588962 


0.589213 


0.589857 


0.601416 


-0.50 


0.691896 


0.693731 


0.702316 


0.734997 
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